Ructuation Statistics and Diffusive Properties of the 2D Triangular 
Lorentz Gas in the Finite-Horizon Regime 



A. Hosseinizadeh" J.F. Laprise'', H. Kroger" '^, G. Melkonyan" '', R. Zomorrodi" 

" Physics Department, Laval University, Quebec, QC, Gl V OAS, Canada ^ 
* Unite de recherclie en sante des populations, CHA universitaire de Quebec, Quebec, QC, GIS 4L8, Canada 
Franlifurt Institute of Advanced Studies, Goettte Universitat Frankfurt, 60438 Frankfurt am Main, Germany 



Abstract 

We study chaotic behavior and diffusion in the 2D periodic Lorentz gas in the finite-horizon regime. The 
dynamical observable which we consider is the length of single particle's trajectories, which moves in a 
triangular array of rigid disks. To analyze the chaoticity of this system, we construct a matrix of the length 
of trajectories and perform a level spacing analysis of the spectrum of this matrix. We find that a universal 
behavior occurs both in level spacing distribution and spectral rigidity. In order to examine diffusion in 
this classical system, we investigate the variance of the length of trajectories versus number of collisions to 
disks. In the case where there is a finite-horizon, such a variance scales linearly with number of bounces. 
This shows that a normal diffusion exists and a central limit theorem is maintained in this regime. 
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1. Introduction 

The 2D periodic Lorentz gas is a model equivalent to a billiard system, in which the billiard ball moves 
in a network of identical hard disks located on a 2D periodic regular grid lH (see Fig. The billiard ball 
(point-particle) alternates between free motion and collision with the disks. In the literature, Lorentz gas 
has been mostly considered either on a rectangular lattice (equivalent to the Sinai billiard on a torus) or 
on a triangular lattice. Here we consider the latter type. The dynamics of the 2D triangular Lorentz gas 
is basically restricted in two regimes: (i) In the finite-horizon regime (2 < o- < 4/ V3, with cr = L/R), 
free paths between colHsions are bounded, because the scatterers are so dense, such that they block every 
direction of motion (note that R is the radius of disks and L is the central distance of two neighboring disks), 
(ii) In the infinite-horizon regime (cr > 4/ VS) there are trajectories, in which the particle can move without 
colliding. The rectangular Lorentz gas (Sinai billiard) also belongs to this class. In the finite-horizon 
regime, the moving particle is confined to a trapping zone (see Fig. |2]l, but in the infinite-horizon regime, 
it can escape to infinity. There is also a dense-packing regime with cr - 2, where the particle is captured 
between three adjacent disks. 

The distinction between the finite and the infinite regimes, defined above in geometrical terms, is also 
manifested in mathematical and physical laws. For example, in the finite-horizon regime, a central limit 
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Figure 1: Geometry of the 2D Lorentz gas, including a scliema of the point particle trapping zone; shown between three adjacent disks 
in a triangular symmetry. 



theorem (CLT) was shown to be true in the case of periodic Lorentz gas flSl- Chernov and Markarian ||6| 
proved that the distribution of the travel time - to give an example of an observable - obeys the CLT, up until 
the /1-th collision (n oo). Therefore, it converges to a normal distribution, in the case where the number 
of collisions goes to infinity. But when the billiard particle moves with constant speed m, the length of 
trajectory A„ = ut„ displays the same behavior, i.e., the distribution of trajectory lengths P{A„) tends towards 
a Gaussian distribution for n — > oo. Moreover, Bunimovich and Sinai pj] demonstrated that the velocity 
auto-correlation function decays exponentially (exp[-«"]), which implies that normal diffusion occurs and 
is accompanied by a finite diffusion coeflicient. This phenomenon has been confirmed by Machta and 
Zwanzig 101 using a numerical calculation of the diffusion constant. In addition, the system converges to 
Brownian motion. 

In the infinite-horizon regime, however, no CLT has been proven yet. Furthermore, in this case there is 
no convergence to Brownian motion f^. Bunimovich jot) has conjectured that the velocity auto-correlation 
function shows an algebraic decay behavior, which suggests that a standard diffusion coefficient does not 
exist. A number of researchers have carried out numerical experiments, in order to determine the decay be- 
havior of the velocity auto-correlation function. In the case of the rectangular Lorentz gas, such studies have 
been performed by Zacherl et al. 1 10 1 , Dahlqvist and Artuso ifTTIl , Zaslavsky and Edelman 1 12 ] , Armstead et 



al. 01311 . and Courbage et al. Ill4ll . These studies show that the velocity auto-correlation function appears as 

(m m ~ c/t , (1) 

where c is a constant. This indicates that the mean square displacement (the variance in position) behaves in 
the following manner: 

<X')~cflogf, (2) 



which was also proven by Bleher [|15|. The same results were obtained for the triangular Lorentz gas by 
Friedman and Martin, Jr 11611 . as well as by Matsuoka and Martin, Jr ||17il . Another chaotic system, in which 



the mean square displacement displays a logarithmic correction, is the Bunimovich stadium billiard UlSll . 



Balint and Gouezel (l^ showed that the distribution of the length of trajectories tends towards a normal 
distribution, where the limit of infinite number of collisions is concerned. In a numerical simulation (not 
presented here) 1I20I1 . we have shown that in the Bunimovich billiard, the variance of length A obeys 

Var(A) oc nlogn . (3) 

One finds that the scaling factor is «log«, instead of the standard expression y/n (see Ref. fl^l for more 
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Figure 2: The 2D triangular Lorentz gas in the finite-horizon regime. Top panel: Regular distribution of nodes xi , ...xn on the boundary 
of the trapping area. Bottom panel: A sample trajectoiy between two boundaiy nodes (a i and .1:3). 
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details). As consequence, there is however no standard diffusion coefficient and no convergence resulting in 
Brownian motion in this billiard. 

In addition to transport properties such as diffusion, the chaoticity of Lorentz gas has also been studied 
both in classical and quantum chaos In particular, it has been shown that in a quantum 

Sinai billiard, statistical properties of eigenvalues disp lay universality. This, in turn, is compatible with 
GOI^ predictions of random matrix theory (RMT) 1I21I1 . But in this work, we are interested in studying this 
subject from a classical point of view. The strategy which we use here, has been demonstrated by Laprise 



et al. in Refs. 1124(1 and II25I1 . By working on a number of classical chaotic billiards (such as Bunimovich 
stadium billiard, cardioid billiard and an optical stadium billiard), and using methods of quantum chaos, they 
have found that universal behavior also occurs in classical chaos. However, in order to use the customary 
techniques of level spacing analysis in quantum chaos, they have computed the eigenvalues of an action 
(length) matri)(0- As demonstrated in next section, such a matrix is constructed from the length of classical 
trajectories of a point-particle, which moves in the billiard zone with a constant velocity. 

As pointed out above, fully chaotic deterministic systems such as Lorentz gas, may display behavior of a 
random system, like Brownian motion, diffusive character and validity of CLT's. Thus, before discussing the 
Lorentz gas, first we take a look at a stochastic billiard. Then we examine the chaoticity of the 2D triangular 
Lorentz gas in the finite-horizon regime, via level spacing analysis of length matrices. And finally, we will 
investigate the diffusion in this system (as a physical property), in terms of trajectory length. 

2. Stochastic Billiard System 

The collision and diffusion of molecules in gaseous or liquid matter can be described by Brownian 
motion, which can be modelled by an open random billiard. In the limit of free path length (Ax — > 0, Af — » 
and (Aji:)^/Af — > const), the model describes diffusion and is mathematically described by the following 
differential equation, having a probabilistic interpretation 



d,p = D{Vfp . (4) 

Let us consider the simple model of a massive particle, e.g., an atom or a molecule, which from time to time 
collides with other atoms/molecules. After each collision, the particle changes its direction and velocity. 
We assume that the new direction and velocity are randomly distributed. In between collisions, the particle 
moves freely, maintaining direction and velocity. The length of trajectory of the random walker, having 
executed A^i-eb collisions, is given by 

A'reb+l ^,,^ + 1 

A = £ Ax, = J] A, . (5) 

/t=l k=l 

We look at the statistical behavior of A for many trajectories, having in common the same number Meb of 
collisions. For statistical purpose it is meaningful to hold Nieb fixed, which means that we keep fixed the 
average length of trajectory 

<A> = (A?,,b + l)<Ax) (6) 



^Gaussian Orthogonal Ensemble 

"^When the particle moves with a constant velocity, the length of the particle's trajectory is equivalent to the corresponding action. 
Thus, instead of an action matrix we construct a length matrix. 
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and A fluctuates around (A). When we plot a histogram of length P(A), we obtain a Gaussian distribution. 
This result holds provided that we consider macroscopic times, i.e. a large number of collisions. Such 
a behavior can be understood from the CLT: In between collisions, any length segment At = Axt of the 
straight trajectory is a random number drawn from some distribution, say P(x). The total length A then is 
a sum of random numbers. Under general conditions on the distribution the CLT says that the sum of 
random variables becomes a Gaussian in the limit of large A^reb- As a consequence, the length of the path of a 
random walker between any given starting point Xin and any end point Xfi follows a Gaussian distribution, if 
we keep fixed the number A^ieb of collisions, and if A^ieb is sufficiently large. If we construct a length matrix 
Aij , where initial and end points are taken from a given set of points xi, ... then all matrix elements 
follow a Gaussian distribution (same variance and mean). Because A^- is real and symmetric, A represents a 
matrix from a Gaussian orthogonal ensemble (GOE). Random matrix theory predicts that the level spacing 
distribution of the eigenvalues from a GOE matrix follows a Wigner distribution. This is similar to the 
statistical approach to quantum chaos, where the level spacing distributions of energy eigenvalues from a 
quantum Hamiltonian are predicted to behave in accordance with random matrix theory (i.e. like GOE in 
the case of a time-reversal symmetric system). According to the BGS|fl-conjecture B211I . this holds when the 
classical counterpart of the quantum system is fully chaotic. Although the random walker is a purely random 
system and not a deterministic chaotic system, such observed behavior of the classical length of the random 
walker (in analogy to the property of Hamiltonian matrix elements in chaotic quantum systems), led us to 
postulate the following hypothesis; With respect to statistical fluctuations, the length matrix may play the 
same role in classical chaos as the Hamilton matrix in quantum chaos. In the following we investigate this 
hypothesis. 

Here let us mention that the BGS-conjecture does not state that the matrix elements itself of a quantum 
Hamiltonian are distributed like a GOE ensemble. The conjecture rather says only that the statistical fluc- 
tuations of the eigenvalue spacings, obtained from the quantum Hamiltonian, are the same as those from a 
GOE ensemble, giving a Wignerian distribution. In other words, it is possible that the matrix elements of 
the quantum Hamiltonian are distributed quite differently from a Gaussian, but nevertheless its level spacing 
distribution is Wignerian. Such a situation, where the distribution of matrix elements is not GOE but the 
level fluctuation statistics is GOE, occurs in nuclear physics. An exarnple is the distribution of the Hamilto- 
nian matrix elements obtained from nuclear shell model calculations {28'\. In this model, there are vanishing 
Hamiltonian matrix elements. This implies that the number of independent matrix elements is much smaller 
than in a random matrix of the same size. However, Brody et al. ll29[l could show that the 2-body residual 
interaction in the shell model yields matrix elements of random character following a Gaussian distribution. 
In particular, they showed that spectral fluctuation properties from such ensembles with orthogonal sym- 
metry are identical to those from GOE. This implies that GOE is meaningful to predict spectral fluctuation 
properties of nuclei governed by 2-body interactions, though the Hamiltonian does not follow a Gaussian 
distribution. 



3. Chaotic Behavior of the 2D Lorentz Gas from the Length Matrix 

3.1. Level Spacing Fluctuations 

In order to study the chaotic property of the Lorentzgas in the finite-horizon regime, first we construct 
the trajectory length matrix A, like in Refs. |[24| and ^3\. For this purpose we consider an array of three 
adjacent disks , as shown in Fig.|2](due to the periodicity of the system). Then we choose a set of boundary 
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Figure 3: Finite-horizon regime (cr = 2.25). Distribution of the length matrix elements for 5000 trajectories. 



nodes {xi,X2, x^}, that are regularly distributed on the border of the trapping region between three disks. 
Because of the symmetry, the boundary nodes are located on one-sixth of the border on the convex section. 

In the next step, we find the trajectory of the point-particle, as it passes between the nodes Xi and xj 
(Fig.|2|. A typical trajectory is constructed by shooting the point-particle from the boundary point x, under 
an initial angle between and n. It moves within the trapping area and following a zig-zag path through 
colliding with the boundary, it finally arrives at the node xj (see Fig.|2l bottom panel). In this way, we are 
able to build a time-reversal invariant, real and symmetric matrix A = [A.jj]f^^f^, in which Ajj is the length of 
trajectory connecting x, and xj. In fact, depending on the initial shooting angle, numerous trajectories can 
be constructed for a specific number of rebounds. According to this, we find an ensemble of length matrices 
for a particular set of boundary nodes. In practice, and to get a better statistics, we perform a level-spacing 
analysis (see below) on each matrix and do a superposition on the unfolded spectra. 

In a given trajectory, indeed, the consecutive segments are correlated to each other However, two typical 
segments with a sufficiently large separation in terms of macroscopic travel-time are statistically indepen- 
dent. Therefore, such segments make an ensemble of random numbers, meaning that the point-particle 
follows a random walk. Through such a random walk, as pointed out before, the length of a trajectory is 
simply given by the summation of free path lengths between successive collisions (i.e., A = 2^ A^ ). In real- 
ity, the so-called step lengths A^ are random numbers from a Gaussian probability distribution function js^l- 
In the case of Lorentz gas we have also obtained the distribution of length matrix elements P(A) which is 
near a Gaussian distribution (see Fig.|3). This means that the matrix A (after subtracting the mean and suit- 
able overall normalisation) approaches a random matrix with GOE statistics 131]. This result is generally 
consistent with the CLT 13]. 

In order to obtain the nearest-neighbor level spacing (NNS) distribution, one has to compute the spectrum 
of the matrix A and unfold the spectrum , 33 1 . This is a customary technique in quantum chaos where the 
number of energy levels tends to infinity. But, here the dimension of the length matrix is finite (it depends 
on the number of boundary points), and therefore, one extracts a limited number of eigenvalues. However, 
in quantum chaos it has been shown that a finite number of levels is also sufficient for obtaining a reliable 
level spacing statistics (see Ref. Issll ). Bearing this in mind, we use a Gaussian broadening Il34ll for the 
purpose of spectral unfolding of the length matrix. Then we c omp ute the NNS distribution P(s), as well as 
the Dyson-Mehta rigidity A3(L) from the unfolded spectrum OSll (see the appendix for definitions). In the 
case of finite-horizon regime, as an example we choose cr - 2.25 and consider a sequence of 40 boundary 
nodes, as well as Nreb - 15 collisions with the hard disks. Fig. |4] shows the NNS distribution P{s) and 
the spectral rigidity A^iL). As one observes, the distribution P{s) and the spectral rigidity Ai{L) follow the 
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GOE statistics, which is predicted by RMT l2 1113211 . Similar results also hold for other values of cr in the 
finite regime. This GOE property establishes universal behavior in the limit of numerous bounces for the 
2D Lorentz gas with a finite-horizon. It is similar to the behavior of energy level spacing distributions in 
quantum chaos, in accordance with RMT. 



3.2. Evaluation of Numerical Errors 
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Figure 5: Top panel: relative error as a function of number of rebounds in finite-horizon regime of the 2D Lorentz gas (for cr - 
as a fully chaotic system. Bottom panel: the same function for a circular billiard as a non-chaotic system. 



2.25) 



In order to estimate numerical errors in the computation of a trajectory, we also compute the time- 
reversed trajectories and compare the starting point of each trajectory with the endpoint of its corresponding 
reversed trajectory. The behavior of such an error in terms of number of rebounds is shown in Fig. |5](top 
panel), with the following observation: 



log<err) = 0.397M-eb - 16 , (M-eb < 40) . 



(7) 



From this relation one observes that the relative error {err) increases exponentially with Meb, that is (err) ~ 
exp[0.397A^reb] (as it relates to the largest positive Lyapounov exponent of the system). Figure|5]displays that 
log(err) saturates around zero for A^reb ~ 40, i.e., there is a hundred percent of numerical error at this point. 
Furthermore, such an exponential form is not observed for integrable billiards. We have examined this for a 
number of integrable systems, and they display a different behavior For example, see Fig. |5](bottom panel) 
for a circular billiard 113 611 . Hence, it can be considered as another fingerprint of chaoticity in chaotic billiards 
such as Lorentz gas. In our calculations we have considered trajectories with A^reb ^ 20, which allows us to 
retain 8 digits of numerical precision for the spectral fluctuation analyses. 
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4. DiiFusion coefficient via tfie iengtfi of trajectory 



Now, let us discuss the diffusion coefficient in terms of the length of trajectory, within a trapping zone 
between three neighboring rigid scatterers (Fig. |2). We have chosen such system, because (i) it is a fully 
chaotic system, having convex repellers/scatterers, (ii) it represents an example showing universality in the 
length matrix A, i.e., Wignerian form in the level spacing distribution P{s), and near a Gaussian behavior in 
the distribution of the matrix elements P{A), and (iii) its diffusion coefficient has been widely studied in the 
literature (for review see [j?! [s?! and references therein). 




< 0.2 



Figure 6: Finite-horizon regime. Top: Variance Var(A) versus number of rebounds for a variety of cr values 
(2.01,2.05,2.10,2.15,2.20,2.25,2.30). Bottom: Diffusion coefficient versus number of collisions for cr = 2.10, 2.25, 2.30. 



As pointed out previously, for the Lorentz gas model in the finite-horizon regime, a number of properties 
have been proven analytically: A central limit theorem holds, the system has the character of Brownian 
motion (in the limit of many bounces), and a diffusion coefficient exists. The diffusive character is manifested 
by a linear relation between the travel time of the point-particle and the variance of position. The diffusion 
coefficient D (in 2 dimensions) is given by the Einstein relation 



At 



(8) 



where t denotes the time of travel and aP denotes the variance of the position (see Ref. Ill5ll ). Now we will 
show that the diffusive character is also apparent via the length matrix which displays universal behavior In 
particular, we claim that there is a linear relation between the time of travel and the variance of the length 
of trajectories with n collisions. That is, AA^ oc t, where travel time t is related to mean trajectory length 
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of n collisions (A„) and velocity u via t 
variable A by 



(A„) /u. We define the diffusion coefficient with respect to the 
/'ar(A) 



(9) 



(note Var(A) = AA^). In numerical simulations, this system was found to display universal behavior man- 
ifested in the spectrum of the length matrix A, as well as in a Gaussian behavior of the distribution P(A) 
(Fig. |3]l. The results for variance of length of trajectories versus A^ieb are shown in Fig. |6] (top panel). We 
performed a statistical average over a number of A^stat = 5000 initial conditions (position and direction). One 
observes that Var(A) scales linearly in A^reb for a variety of the geometry parameter cr in the finite-horizon, 
ranging from 2.01 to 2.30. This linearity is in a general agreement with the results of Borgonovi et al. llssll . 
An estimate of diffusion coefficient (Eq.|9]l is also shown in Fig. |6] (bottom panel) for three different values 
of cr. From this figure one observes that ©a's are relatively constant. 





Figure 7: Top panel: The moving zone between adjacent disks in the case of the localized regime of the 2D Lorentz gas. The point 
particle is trapped in this zone. Bottom panel: Diffusion constant of the Lorentz gas in the localized regime with o" = 2.0. 

In the limit of the dense-packing (or the so-called localized) regime (with cr = 2), on the other hand, we 
found Da ~ (Fig.]?). This was expected, because the particle is caught between connected disks with no 
opportunity to escape beyond the trapping area. 

Finally, when going from the finite (2 < cr < 4/ V3) to the infinite-horizon regime (cr > 4/ V3), it has 
been predicted that the periodic Lorentz gas undergoes a change in behavior [10, .15, 3X,i390. In the infinite- 
horizon regime, the standard diffusion coefficient D is known to diverge iH. Normal diffusion goes over to 
anomalous diffusion, where the variance of position scales like the number of rebounds with a logarithmic 
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Figure 8: Distribution of lengtli matrix elements P(A) in the infinite-liorizon regime with o" = 2.5. Number of rebounds and number of 
boundary nodes are A'reb = 15 and N = 40, respectively. 

correction (as mentioned in Eq. ([3)) ifisll 

aP oc n login) . (10) 

We have done a numerical study of the triangular Lorentz gas in the infinite horizon regime, searching for 
the scaling law of the form AA^ oc « log(«). However, our data were too noisy, being due to long trajectories 
in between bounces with the disks. The reason of such trajectories is that in this regime, the point-particle 
finds some corridors and travels long distances without colliding to hard disks. Because of this, even the 
length distribution P(A) differs from a Gaussian, as shown in Fig.|8] With presently available computational 
resources we were not able to reduce numerical errors, and to obtain statistically significant results to test 
such a logarithmic behavior This can be considered as a future research work. 

5. Summary and Remarks 

In this work we have studied the chaotic behavior and diffusive properties of the classical 2D triangular 
Lorentz gas, using the length of trajectories. Based on RMT and through level spacing analysis of length 
matrices, we have succeeded in finding the following results: 

(i) The NNS distribution P{s) and the spectral rigidity A3(L) (which are obtained from eigenvalues of the 
length matrix A) are in conformity with a Wigner distribution. This is consistent with the GOE prediction 
of random matrices. 

(ii) The distribution P( A) of length matrix elements resembles a Gaussian distribution. This is in agree- 
ment with the random walk and Brownian motion, as direct consequences of the CEP. In the case of infinite- 
horizon, the distribution P( A) deviates considerably from the normal distribution. This is possibly related to 
the long free paths of the point-particle. 

(iii) The numerical errors (in the finite-horizon regime) exponentially evolve as a function of the number 
of rebounds. Such a behavior does not exist in integrable billiards, and therefore it can be considered as 
another signature of chaos in the system in question. 

(iv) The variance of the length of trajectories linearly varies in terms of the number of collisions, which 
yields a normal diffusion coefficient when finite-horizon is involved. In the localized regime, there is no 
diffusion because the point-particle cannot move between rigid disks. 



I ^Data I 
Normal fit 




4 6.8 10 12 

A 



11 



Finally, we should notify that the works establishing universality in semi-classical systems are usually 
based on classical periodic orbits (for example, see Ref. H). But, here we showed that universality can 
also be observed using non-periodic bouncing trajectories. 

Acknowledgments. Helmut Kroger has been supported by NSERC Canada. 



6. Appendix: Spectral Statistical Analysis 

Consider a quantum system with the spectrum e, (; = 1,2, ...A^) and the level spacing i, - e,+i - e, 
between subsequent levels. In general, the spectrum distribution N{e) of a Hermitian matrix fluctuates with 
its eigenvalues. One can define a mean level distribution N{e), and therefore, the primordial level density 
A^(e) can be sphtted into an average part, as well as a fluctuating part: 

Nie)^Nie) + Nnucie). (11) 

The average level distribution is given by 

N{e)^ [ de'p(e'), (12) 



where p(e) is the mean level spacing density function. 

In fact, N{e) is the dominant system dependent part, and Na^de) exhibits the universal part of level spac- 
ing distribution. One performs such a separation because the interesting part of the level spacing distribution 
is the universal fluctuating part, which takes a small contribution of the distribution. In a fully chaotic system. 



one expects that this part is independent of the particular system, and exhibits universal fluctuations 1133 1. 

The main goal here is to study the sub-dominant fluctuating part of the spectrum by RMT. For this 
purpose, first one needs a constant mean level spacing s of eigenvalues. The technique that people use in 
quantum chaos is to renormalize the eigenvalues, in such a way that their average separation density p(e) 
becomes unity. After this so-called unfolding procedure, the leading smooth part N{e) gets a unit spacing on 



statistical average. The unfolded levels are characterized by 113411 



{w,) = {A?(e,)) , /=1,2,...A?, (13) 



where w,- are dimensionless levels with the separation 5, = - w,, and with the density p(w) = 1. For the 
unfolding process one can use a Gaussian broadening method, where the average level density function p(e) 
is given by I34i] 

{e-e,f 

y pxn - 

•V2^- 



2cr2 



(14) 



From the unfolded spectrum, one computes the NNS distribution P(s), which denotes the histogram of 
finding a separation s between neighboring levels. For a chaotic system it is a Gaussian, i.e., Pwigner('S) - 



nsl2 exp[-7ri74] , ^ > 14111 



When the rescaled sequence of levels is available (i.e. {w, | / = 1,2, ...,A^)), then one can also evaluate 
the least square deviation of the spectral step function N{<x>), from the best straight line (i.e., au + b) fitting 
it in an interval [a, a + L\of the length L. It is called the Dyson-Mehta rigidity (or spectral rigidity), given 
byil 

Ai(L; Q-) = - Mina^b I [N{aj) -au- bf dw . (15) 

^ Ja 
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By replacing the functional form of N{a>) in Eq.lfTSll. one obtains a formula as 11351] 

-i2 



A3(L; a) 



16' I? 



3n_ 
21? 



(16) 



3_ 
1? 



+j ^(«-2/+ l)w,- 



with cji - (Di - (ff + |), which means one takes the center of the interval as origin. Note that other measures 
of the spectral rigidity, i.e. Ai and A2, have been also introduced by Dyson and Mehta ll32ll . but A3 has 
proven to be more useful in RMT. 
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